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Abstract 

We derive integrable discrete systems which are contiguity relations of two equations in 
the Painleve-Gambier classification depending on some parameter. These studies extend 
earlier work where the contiguity relations for the six transcendental Painleve equations 
were obtained. In the case of the Gambler equation we give the contiguity relations for both 
the continuous and the discrete system. 

t Permanent address: CRM, Universite de Montreal, CP. 6128, Succ. Centre-ville, Mtl, II3C 3J7, Canada 



1. Introduction 



The relation of discrete Painleve equations to their, better known, continuous analogues has 
by now been firmly established. As a matter of fact one of the very first methods proposed 
to the derivation of discrete Painleve equations (d-P) has been the one based on the auto- 
Backlund/Schlesinger transformations of continuous Painleve equations [1]. (This link goes 
even further back in time since one of the first discrete Painleve equation obtained, even 
before the notion of a d-P was explicited, was the one derived by Jimbo and Miwa [2] as 
a contiguity relation of the solutions of Pn). The general method for this derivation has 
been amply explained in previous publications [3]. One uses the auto-Backlund/Schlesinger 
transformations in order to derive relations for the solutions of a continuous Painleve equa- 
tions for contiguous values of the one of its parameters and one gets precisely a d-P. The 
discrete equation related to Pn has been first obtained, as explained above, by Jimbo and 
Miwa [2] and rediscovered in [1]. The d-P's related to Pm were first given in [1] and have 
been recently reexamined in detail in [4]. The discrete equation related to Piv was first 
obtained in [1] and was again studied in [5]. Some d-P's related to Py were first given in 
[1] but the relation of Py to asymmetric d-Pn is under preparation [6]. Finally the discrete 
systems obtained from the auto-Backlund of Pyi were presented in full detail in [7] . Thus 
the discrete, difference, Painleve equations can be interpreted as the contiguity relations of 
the continuous ones. (This would tend to cast a doubt on the fundamental character of the 
d-P's. While this is justified as far as the difference d-P's are concerned there exists a second 
type of d-P's which have no relation to continuous Painleve equations whatsoever. These 
are the multiplicative equations, g-P's [8] which are purely discrete objects and thus of the 
most fundamental nature). One of the consequences of this relation between continuous 
and discrete Painleve equations is that one can obtain the Lax pairs of the latter from the 
mere existence of the Lax pairs and Schlesinger transformations of the former. The only 
complication comes from the fact that given a d-P one has to identify the procedure for its 
derivation based on the Schlesinger transformation of some continuous Painleve equation. 
This can be done in the frame of the geometric description of Painleve equations which was 
introduced in [9] and dubbed the "Grand Scheme" . This description based, on the geometry 
of the affine Weyl group weight spaces, has turned out to be a most valuable tool for the 
classification of the d-P's, which minimised their uncontrollable proliferation. 

In the approach we have described above the emphasis was on the discrete analogues of the 
six transcendental equations of Painleve. As a matter of fact the remaining integrable ODE's 
of the Painleve-Gambier classification have received little attention. Our current knowledge 
is essentially limited to the discrete analogues of the 24 basic equations of Gambler [10]. 
In particular no results exist on the contiguity relations other than for the six equations 
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of Painleve. In this paper we intend to focus on just this question and present a first, 
exploratory, study of some discrete equations related to continuous systems which are not 
one of the six Painleve equations. In particular we shall obtain the mappings related to the 
equations 35 and 27 of the Painleve- Gambler classification. Equation 27 is also known as 
the Gambler equation and has been studied in detail in [11], where its discrete equivalent 
has been obtained. Since, as was shown in [12], both the continuous and discrete Gambler 
equations possess Schlesinger transformations, one can obtain the contiguity relations not 
only of the continuous but also of the discrete system. This result is presented in the last 
part of the present paper. 



2. Discrete systems related to P35 



The first system we are going to examine is the equation #35 in the Painleve-Gambier 
classification [13]: 

w" = —to - \ ^w- -q )w' + -W^ - —qw^ + { 4q' + r + -q^ ] w + 2qr - 3r' 

Sw \S S w J S 6 \ ^ ) ^ 

(2.1) 

with r = -^/3 - 2(5' + and q one particular solution of Pn, 

v" = 2v^ + zf; + a, (2.2) 

where z is the independent variable and a is a constant. The way to integrate this equation 
is based on the following Miura transformation. We start with the general solution of Pn w, 
and we construct the quantity 

w = ^ ~^ + V + g (2.3) 

V — q 

which solves (2.1). Moreover, starting from the general solution of (2.1) to, the quantity 

j^w-i - z) (2.4) 

solves Pn- Equations (2.3) and (2.4) define a Miura transformation relating Pn and P35. 
Another such Miura relation between Pn and some other equation of the Painleve-Gambier 
list is already known, that relating Pn and P34. It has been studied in detail in [14] in both 
the continuous and the discrete case. 

Equation P35 has a parameter, a, through its relation to Pn- Let us write for the solution 
of (2.1) corresponding to the parameter a and and q^ for the solutions of Pn entering in 
(2.1), (2.3) and (2.4). Next we consider the auto-Backlund transformations of Pn [15]: 
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corresponding to solutions where the parameter a varies by integer quantities. Thus, fol- 
lowing the general theory presented in [1] one can introduce a discrete system associated 
to these transformations as a mapping which is the contiguity relation for the solutions of 
Pii- This equation was in fact first obtained in [2] and is now known (at least to the present 
authors) under the name of alternate d-Pj. It has the form 

a - 1/2 Q + 1/2 2 



+ — — = -2voc - z. (2.6) 

Va-l+Va Va+l+Va 

This is the equation satisfied by all solutions of Pn, and in particular by q, when we vary 
the parameter of Pn using (2.5). 

Using the Miura (2.3) we can obtain the solutions of P35, Wa+i and Wa-i, corresponding to 
the parameter values ail: 



Wa+1 = h Va+1 + Qa+l, (2.7a) 

Va+1 — Qa+l 

v' — n' 

Wa-1 = h Va-1 + Qa-l- {2.1b) 

Va-1 - Qa-l 

Using (2.2), (2.3), (2.5) and (2.7), one can write an equation relating Wa+i, Wa and Wa-i 
alone where neither v nor its derivatives appear. Interpreting a as a discrete independent 
variable, we obtain thus a second-order discrete equation for w. This mapping is however 
quadratic in Wa+i and Wa-i and thus (as was also argued in [16]) cannot be integrable. 
Our argument is based on the fact that the evolution of such a mapping leads in general to 
an exponential number of images, and preimages, of the initial point. The non-integrability 
of this non-singlevalued system is not in contradiction to the fact that we can obtain one 
solution of this mapping, namely the one furnished by the evolution of the system before 
elimination. This solution is the only one that we know how to describe, while a system 
with exponentially increasing number of branches eludes a full description. 

In order to proceed, instead of working with a second-order mapping, we choose to derive 
a system of two coupled equations. We first make use of the discrete symmetry of the 
solutions of Pii: if Va is a solution of Pn with parameter a, then —Va is a solution of Pn 
with parameter —a. We then introduce the quantity 

Ua = ^^^^^ -v^-qa (2.8) 

Va - Qa 

which is obtained from the expression (2.3) for Wa after implementing the reflections Va — > 
—Va and Qa —Qa- The quantity Ua is thus the solution of a slightly modified P35 equation 
(with —Qa instead of Qa) that we shall denote P35. Note that, from (2.3) and (2.8) we have: 

Wa-Ua^2{Va + qa)- (2.9) 
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We start by writing an equation between Wa-i, and Ua only, using (2.2), (2.3), (2.5), 
(2.7) and (2.9). Similarly we also obtain an equation relating Ua+i, and w^- We thus 
find the second-order discrete system: 

= + ■ 2^—7 r . , (2.10a) 

Wa^Uoc^ ■ r H , (2.106) 

Ua{qa + qa-l) Wa-1 - 

where is a particular solution of (2.6). 

This system introciuces the contiguity relations for equation P35 (2.1) and P35. It can be 
integrated with the help of alternate d-Pi, equation (2.6). We first write v'^ in terms of Va+i 
and Va and q'^ in terms of Qa+i and using the Schlesinger transformation (2.5). Using 
these expressions in the definition (2.3), we obtain the following: 



a + 1/2 

Va - qa 



qa + qa+1 Va + f a+l 



(2.11) 



The quantity is obtained in terms of Va and Wa using (2.9). Thus, once Va is known, we 
can construct Wa and in a straigthforward way. 

At this point we can remark that while (2.10) is integrable by construction, the condition for 
integrability, namely that g is a specific function, may appear somewhat contrived. In what 
follows we shall show that the fact that q satisfies alternate d-Pi can be obtained precisely 
through the singularity confinement criterion. We start by writing (2.10) and (2.6) in a 
more convenient way in order to implement the singularity analysis.: 

Ua+l^Wa + -, rj ■ r r— r , (2.12a) 

Wa^Ua 7 ■ tH , (2.126) 

Ua{qa + qa-l) Wa-l - Ua 

where q and z are two functions to be determined. 

There are three possible sources of singularities for (2.12): Ua — 0, Wa — and the case 
when the denominator in the rhs of equation (2.12a) happens to be for some value of 
a. We first examine the singularity Ua — 0. Let the initial condition Wa-i be free and 
introduce the small parameter e: — e. At the lowest order in e we find, Wq, = 1/e + . . ., 
Ua+i — 1/e + . . ., Wa+1 — 0{e). Then Ua+2 must take the indefinite form 0/0 in order to 
contain the information on Wq-i when e ^ 0. Implementing this condition, we find that q^ 
must indeed satisfy 

^ + - -2qJ - 7, (2.13) 



9a- 1 + qa qa+1 + qa 
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for some constant 7 and Za must be linear in a. (The two other singularities are automati- 
cally confined for all functions qa. and Za.) Thus the singularity confinement criterion gives 
for Qct and z^ precisely the constraints expected from the construction of (2.12) from P35. 

3. A MAPPING RELATED TO THE CONTINUOUS GaMBIER EQUATION 



y' = -2/^ + C, 
x' — ax^ + nxy + a, 



The second discrete system we are going to derive is the contiguity relation of the solutions 
of the Gambler equation. The continuous Gambler system can be written as a cascade of 
two Riccati equations: 

(3.1a) 
(3.16) 

where a, c and a are free functions of the independent variable z and n is an integer. By a 
change of variable, the function a can be scaled to or 1 but, for the purpose of constructing 
Schlesinger transformations, it is more convenient to keep it as a free function. Moreover, if 
we require that equation (3.1) possess the Painleve property, the functions will be subject 
to one more constraint. This constraint depends on the value of the integer parameter n. 

Let us write Xn, an and an, the quantities appearing in (3.1) corresponding to the coupling 
parameter n in (3.1b). Then, we can introduce the following Schlesinger transformation 
[12,17]. It relates the solutions of (3.1) to that of an equation of the same form but with a 
coupling parameter n — n + 2. This transformation is given by: 

-1 



Xn+2 — 



'n+2 



O'n+2 

n + 1 

n 



y + 



+ 



a>nXr 



n + 2 ttn, n + 1 



C"n+2 
O'n+2 ■ 



n 



2 ttr 



n 



1 



Cn + 



ttnCTr. 



n 



3 a 



12 



0"n+2 V n+1 

The inverse of this transformation is given by 

1 — n 



n 



2 Or 



n 



2)2 a: 



Xn-2 = 



''n-2 



Ct.n-2 

n 



(Tr 



Xnin - 1) 



+ y + 



n Ur 



O'n-2 



n 



2 Or 



0"n-2 = 



1 — n 

Q-n-2 



Cn + 



+ 



+ 



n — 3 



J2 



n 



n Or, 



{n - 2f 



(3.2o) 
(3.26) 
(3.2c) 

(3.3a) 
(3.36) 
(3.3c) 



We can obtain a second-order mapping for x eliminating y between equations (3.2a) and 
(3.3a). 



Xn-2- 



O-n-2 



ttr 



n 



1 1 

+ 



(^n+2 



n + 1 (T„n — 2 a„n-|-2 Xn.n — 1 Xn+2 'n+ 1 







(3.4) 
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There is another way to present equation (3.4) if we write it as a system, involving expHcitly 
y and using the fact that it does not depend on n. The two equations (3.2b) and (3.2c) 
are equivalent to the upshifts (by 2 units of n) of (3.3b) and (3.3c). Equation (3.2b) is the 
compatibility condition between (3.2a) and the upshift of (3.3a). 

We obtain thus 

yn+2 = Vn (3.5a) 

^ (3.56) 



Xn-2 = 



Or 



<Jr,. I 

+ 2/n + 



Xnin — 1) 2 — n CTr^ 

The evolution of a and a is obtained through equations (3.2b)-(3.2c). System (3.5) is just 
a particular case of the Gambler mapping introduced in [18]. It is remarkable that the 
contiguity relation of the solutions of the Gambler equations with respect to the one integer 
parameter, namely n, is again a (discrete) Gambler equation. However the latter does not 
have the full freedom of the Gambler mapping (in perfect parallel to the relation between 
continous Painleve equations and d-P's.) 

4. A MAPPING RELATED TO THE DISCRETE GAMBIER EQUATION 

In this section we shall concentrate on the generic discrete version of the Gambler equation. 
As was shown in [12] it can be written as a system of two discrete Riccati's in cascade: 

Vn+i = - — — -j-' (4- la) 

XniUn "^n) ~l~ QniUn ^n) -. 7 \ 

Xn+1 — • K'^-^'^) 

XnUn 

(Other choices for the Gambler mapping do exist but equation (4.1) is the most convenient 
for the derivation of the Schlesinger transformation.) The mapping was analysed in [17] 
using the singularity confinement method. The only singularity that plays a role is that 
induced by a special value of y, which, given the form of (4.1), is y — 0. The way for this 
singularity to be confined is to require that the value of x go through an indeterminate value 
0/0 after a certain number of steps, say N . This means in particular that x must be zero, 
while s appearing in (4.1b) must be equal to the quantity ipjsi (introduced in [12]) which 
is the A^th iterate oi y — Q under (4.1a), N times downshifted i.e. V'o.n — 0, V'l.n — Cn-ii 
'^2n — '^""^"'"^"^ , etc. As was shown in [12], at the continuous limit, the number N of 
steps, necessary for confinement, goes over to the parameter n appearing in the continuous 
Gambler equation (3.1). 

The discrete version of the Schlesinger transformation (3.2) is a transformation which relates 
the solution of a mapping confining in N steps to that of a mapping confining in A^ + 2 steps. 
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The complete analysis for the derivation of the discrete Schlesinger can be found in [12]. 
Here we shall just state the result: 

Xn+1 = <Pn 1 7 ry— , (4.2) 

where x is the solution requiring + 2 steps for confinement and x stands for the solution 
necessitating just steps. The quantity is given by 

^ ^ (V^JV,n-l + l)(^Ar,n-lV^Ar,n-2gn-l + i>N,n-2rn-\ " 

,nQn + rn) + 

With this proper choice of we find that x satisfies an equation: 

iUn — ^n) + QniVn — '4^N+2,n) 
Xn+1 — Z ) V^-^l 

Xnt/n 

where f and q are functions completely determined in terms of g, r and c. In a similar way 
we can derive the Schlesinger transformation invoving a solution which takes N — 2 steps to 
confine. We find: 

XniUn ^n— 1 

Xn — '^n 7 \ , , / ; x") V^-^) 

XniVn - Cn-l) + KijJn - WN,n) 



where 



,n+l 



(4.60) 



Cn-l + C„ - 1pN,n+l{Cn-l + I) 
^ ^ ttn(l - Cn) + Kjl + Cn - 2-i/'jv,n+l) g x 

Wn(l - Cn) + fl-nCl + Cn - 2V'jv,n+l) ' 

«n = l-rn + ^^%^. (4.6d) 
/^n-l 

We can now obtain the contiguity relation of the solutions of the Gambler mapping by 
eliminating y between (4.2) and (4.5). 

On the other hand, just as in the case of the continuous Gambler system there exists a 
simpler way to write this mapping. The value of y appearing in (4.2) or (4.5) is independent 
of N. Thus a far simpler way to present the mapping is for instance as 

Vn^Vn, (4.7a) 

^ XniVn - Cn-l) + 5'n(|/n " ^JV,n) ,^ 

''Xn{yn-Cn-l)+K{yn-'lpN,n)' 

Again we remark that (4.7) is a mapping of Gambler type although of a very particular one. 
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5. Conclusion 



In this paper we have studied the contiguity relations of the solutions of two equations in the 
Painleve- Gambler classification which contain some parameter. This kind of investigation 
was previously limited to the six transcendental equations of Painleve. The main result 
of those studies was to show that the discrete Painleve equations are just the contiguity 
relations of the continuous Painleve equations. Here we have extended this approach to 
equations which, while not being one of the six Painleve equations are integrable and, in 
fact, belong to the same classification. In the case of the Gambler equation we have used the 
Schlesinger transformation which is associated to the integer parameter n appearing in the 
equation. Our result in the case of the continuous Gambier equation is that the contiguity 
relation of the solutions satisfies a mapping which assumes the form of a discrete analogue 
of the Gambier system, in a nice analogy to the situation for the continuous and discrete 
Painleve equations. In parallel to our approach for d-P's the discrete Gambier system has 
also been analysed. We have shown that using the Schlesinger transformation that also exist 
in this case one can derive the contiguity relation of the solutions which turns out to be in 
the form of (a special case of) a Gambier mapping. We expect our approach to be applicable 
to other equations of the Painleve- Gambier classification which contain parameters and for 
which Miura/Schlesinger transformations exist. Now that the question of classification of 
discrete Painleve equations is in the process of being settled [19,20] it is interesting to analyse 
the remaining equations in the classification and work out their implications in the discrete 
case. 
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